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Abstract. Theoretical models have had difficulty matching the obse rved number 
densi ty of sub-millimeter galaxies (SMGs), causing some authors (e.g., iBaugh et al.l 
2005) to suggest that SMGs provide evidence for a top-heavy initial mass function 
(IMF). To test this claim, we have, for the first time, combined high-resolution 3-D 
hydrodynamic simulations of isolated and merging massive, gas-rich galaxies, radia- 
tive transfer, and a semi-empirical merger rate model to predict the number density of 
SMGs. Our model can reproduce the observed SMG number density even when using 
a standard (Kroupa) IMF. The agreement is due to a combination of relatively long sub- 
mm duty cycles for mergers (a few times 10 8 years for our most massive models), which 
owe to our combination of high-resolution 3-D hydrodynamic simulations and dust ra- 
diative transfer; sufficient number densities of massive, gas-rich me rgers; and the de- 
creas e in sub-mm counts observed by recent deep/wide surveys (e.g. jAustermann et al.l 
1201 Oh relative to previous surveys. Our results suggest that the observed SMG number 
counts do not provide evidence for a top-heavy IMF at high redshift. 



1. Introduction 

Sub-millimete r galaxies (SMGs) are a population of extremely luminous (Lboi ~ 10 12 - 
10 13 L Q ; e.g.. iKovacs et all [20061) high-redshift gal axies that emit a lmost all of their 



bolometric luminosity in the IR (for a review see iBlain et al. 2002). SMGs are pri 



marily powered by sta r formation rather than AGN (C. Hay ward et al. 2010a, in prep.; 
Alexander et alj|2005l) . with inferred star formation rates (SFRs) of order 10 2 - 10 3 M Q 



yr . Some models suggest that the stellar initial mass function (IMF) may be top- 
heavy in such high-density environments (e.g., Mark Krumholz's proceeding in this 
volume), so SMGs are an excellent place to look for IMF variation. 

The sub-mm flux of a galaxy depends on the amount of energy absorbed by dust 
and the mass, geometry, and composition of the dust. A top-heavy IMF means that more 
massive stars are formed per unit SFR, resulting in greater luminosity emitted and dust 
mass produced per unit SFR. The combination of increased luminosity heating the dust 
and increased dust mass leads to greater sub-mm flux. Therefore the sub-mm flux of a 
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galaxy - and thus the number density of sub-mm sources - predicted by a theoretical 
model will depend sensitively on the assumed IMF. 

To date, models have had difficulty reproducing the observed sub-mm number 
counts usi ng a standard IMF . The number counts predicted by the semi-analytic model 
(SAM) ofBaughetaD|2005) were less than the observed counts by factors of > 20. Of 
the model variations Baugh et al. tested, only by assuming a flat IMF (dN /d log M = 
constant) in starbursts were they able to reproduced the observed counts. This increased 
the sub-mm counts from starbursts by a factor of ~ 10 3 . However, the large parameter 
space of SAMs often yields multiple qua litatively distinct solutions that satisfy all ob- 
servational constraints (e.g. jLu et al.ll201fjh . leaving open the possibility that the SAM 
of Baugh et al. does not requir e a top-heavy IMF to match the observed counts. 

As an alternate approach, iDave et all (1201 fj) presented a large cosmological sim- 
ulation for which they used the observed SMG number density to calculate how many 
of their simulated galaxies should be SMGs, assuming that the SMGs were the most 
highly star-forming galaxies in the simulation. The galaxies they identified as SMGs 
had SFRs < 3 times those inferred from the IR luminosities of SMGs. One possible 
explanation for the discrepancy is that the IMF in SMGs is top-heavy; if so, the star 
formation rates inferred assuming a standard IMF would be higher than the true SFRs 
of the SMGs. (Note that the required change to the IMF is much milder than that re- 
quired by Baugh et al.) However, there are multiple other potential explanations for 
the discrepancy (see Dave et al. for discussion of these). Furthermore, it is not clear 
if the resolution of the Dave et al. simulation is sufficient to fully resolve the SFR en- 
hancement caused by tidal torques in galaxy mergers. Finally, our models show that the 
simple rank ordering between SFR and sub-mm flux assumed by Dave et al. does not 
hold (C. Hayward et al. 2010b, in prep.), so the evidence provided by their work for a 
top-heavy IMF is indirect at best. 

The discrepancy between modeled and observed SMG number densities provides 
tantalizing but indirect evidence for a top-heavy IMF. However, the complexity of the 
modeling leaves open the possibility that something other than a top-heavy IMF can 
explain the discrepancy. In this work we present a novel method to predict the sub-mm 
number counts: we combine high-resolution hydrodynamic simulations of isolated and 
merging disk galaxies, 3-D radiative transfer (RT) including a fully self-consistent cal- 
culation of the dust temperatures, and a semi-empirical model for merger rates. Prelim- 
inary results show that our model can reproduce the observed sub-mm number densities 
even when a standard IMF is used. 



2. Methodology 



In iNarayanan et al.l (1201 Pel) we presented a merge r-driven model for SMGs. We per- 
formed high-resolution Gadget-2 (ISpringell 2005) 3-D hydrodynamic simulations of 
merging and isolated massive, gas-rich disk galaxies scaled to z = 3 and used the 
3-D adaptive grid polychromatic Monte Carlo dust RT code Sunrise (IJonsson 2006; 
Jonsson et all 120101) to calculate the spectral energy distributions (SEDs) of the simu- 
lated mergers. The model successfully reproduced the range of sub-m m fluxes, UV-mm 
SEP, and bulge, black hole, and gas masses of SMGs. Furthermore, in lNarayanan et al. 



(2009) we calculated the molecular line emission o f the simulated galaxie s using the 3- 
D non-LTE molecular line RT code TurtleBeach ( Narayanan et al. 120081) . finding that 
our models were consistent with the observed line widths, compact CO spatial extent, 
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and excitation conditions of SMGs. Finally, our models naturally explain the inter- 
section of the SMG and dust-obscured galaxy (DOG) populations dNarayanan et all 
2010bh . While the consistency between our models and many observed properties 
of SMGs is encouraging, we have not yet shown that our models can reproduce the 
observed SMG number counts, a key test of potential SMG models. Calculating the 

number counts predicted by our models is the goal of the work describ ed here. 

He re we recapitulate the key details of the simulations presented in Narayana n et al. 
(2010c) and describe the primary difference between our simulations and theirs. The 
Gadget-2 hydrodynamic simulatio ns include star formation via the volume-density- 
dependent Kennicutt-Schmidt law ( Kennicutt 1998 : S^hmidt|[l95 9) (consistent with ob- 
servations of high-redshift galaxies; e .g., Narayanan et ~ 2010ah . the two-phase ISM 
model of ISpringel & Hernquist (2003, hereafter SH03), and black hole accretion and 
feedback following ISpringel et al.l (I2005h. I n post-processing we use the dust RT code 
Sunrise ( Jonssonll2006l : Jonsson et al.ll2010h to calculate the UV-mm SEDs of the sim- 
ulated galaxies assuming a lKroupal d200lh IMF. Sunrise calculates the temperature of 
each dust grain species assuming thermal equilibrium and self-consistently treats dust 
self-absorption, which is key for an accurate determination of the sub-mm flux owing 
to the high optical depths of the central regions of the me rgers. 



The calculations presented here differ from those of iNaravanan et al.l (I2010cl) pri 



marily in the treatment of the obscuration from cold clouds in the ISM. INaravanan et al. 



(1201 Pel) assumed that the young stars present in the simulations were completely ob- 
scured by their mo lecular gas birth clouds for 10 Myr (/pdr = 1 in the parlance of 
Groves et all 12008). The attenuation and du st emission owing t o thes e birth clouds 
was calculated via the sub-resolution model of G roves et al. I (l2008l). (SeelJonsson et all 
(2010) for details of the Sunrise implementation of the model.) INaravanan et a l. (2010c) 
neglected any other obscuration from the cold phase of the SH03 two-phase ISM model, 
so that stellar radiation only encountered dust in the birth cloud (if the radiation is from 
a star particle with age < 10 Myr) and the diffuse phase of the ISM. Instead, we opt to 
neglect any sub-resolution obscuration from the birth clouds (/pdr = 0) and model the 
high obscuration by ignoring the SH03 multiphase breakdown for the purposes of the 
RT calculation, assuming that the total SPH gas density is distributed uniformly on sub- 
resolution scales (< 300 pc). By eliminating use of the SH03 sub-resolution model we 
can more effectively test how the assumptions of our model affect the predicted number 
counts. The consequences of this choice will be discussed in detail in C. Hayward et 
al. (2010c, in prep.). 

In order to predict SMG number counts from our models, we combine sub-mm 
lightcurves fro m the Gadget+Su nrise simulations with the semi-empirical merger rate 
model of iHopkins et all (l2010al) . The number density of sources (deg -2 ) with flux 
density greater than S a is 



dN(> 5a) 
d£l 
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dN 



dVdtd log Mbdjidf^ 



(M b ,n,f g ,z) 



dV 



r(> S a , M b , [i, f g , z)——(z)d log M b d{idf g dz, 
dlldz 



where the term on the right-hand side first line is the comoving number density of 
mergers of fixed most-massive progenitor baryonic mass Mb, merger mass ratio fi, gas 
fraction at merger f g , and redshift z; t(> S A,Mb,ju,f g ,z) is the amount of time (duty 
cycle) a merger with given Mb, >u, and f g put at redshift z has observed-frame flux 
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density at wavelength A greater than S and dV/dQdz{z) is the comoving volume per 
square degree per unit redshift at re dshift z. We integra te over the redshift range 1-6 
and assume a WMAP-7 cosmology (lLarson et al.ll2010h . 

For each of our simulated galaxies we calculate the duty cycle for various sub- 
mm fluxes. We fit the model duty cycles as a function of progenitor baryonic mass 
and merger mass ratio. We neglect any dependence of the duty cycle on gas fraction 
or redshift, as all progenitor disks are initialized with 80% gas (so that sufficient gas 
remains at merger coalescence) and scaled to z = 3, and the negative ^-correction causes 
the observed sub-mm flux of a galaxy to be effectively independent of its redshift within 
the redshift range with which we are concerned. Duty cycles in hand, w e calculate the 
merger rates using the semi-empirical model of Hopkins et all (l2010alJbl) . The model 
starts with an observed stellar mass function (MF) (IMarchesini et al.ll2009h . assigns gas 
masses from observations, and places the galaxies in (sub)halos via a halo occupation 
distribution approach. Given halo-halo merger rates from A^-body simulations, galaxy- 
galaxy merger rates are calculated by assuming the galaxies merge on a dynamical 
friction timescale. 



3. Results 

Figure Q] shows the AzTEC 1.1 mm number counts from mergers predicted by our 
fiducial model, where we have discarded all snapshots with gas fraction greater than 
60% when calculating the sub-mm duty cycles. The predictions agree well with the 



Austermann et all (12010 ) SHADES data, which we use because to our knowledge this is 



the largest field for which deep ~ 1 mm counts are availabl e. Previo us 850 /mi SCUBA 
counts were biased high by cosmic variance (lAustermann et al.l2010l) : SCUBA-2 should 
remedy this problem. One possible reason we do not have any galaxies with 5yi mm > 
5.5 mJy is because we do not extrapolate beyond our largest simulated merger (a ma- 
jor merger of two disks with baryonic mass 4 x 10 u M o ), and the ~ 5 sources with 
S l.i mm > 5.5 mJy may be even more extreme, rare objects. Alternately, the discrep- 
ancy may owe to our underprediction of the sub-mm flux at merger coalescence for the 
reasons outlined below, as a ~ 10% boost of the flux of our brightest sources is enough 
to account for all but the single brightest SMG observed in the SHADES field. 

In our models isolated disks contribute negligibly to the 'classical' SMG popula- 
tion (S 850 fim > 5 mJy); a 60% gas disk with baryonic mass 4 x 10 1 1 M© has S sso fim ^ 4 
mJy (5 i.i nun ^5 1-5 mJy). Mergers contribute to the SMG population in two ways: 
First, as has been suggested by numerous authors, some SMGs are merger-driven star- 
bursts. In this scenario, tidal torques at merger coalescence induce a strong starburst 
and peak in AGN activity, both which cause the luminosity to sharply increase. As 
a result, the dust luminosity and thus sub-mm flux increases, thoug h this is partially 



mitiga ted by an increase in the dust temperature. This is the scenario iNarayanan et al. 



(2010c) studied in detail. The second contribution is from the merging disks before 
they are strongly interacting. The beam sizes of SCUBA and AzTEC are ~ 30", which 
is ~ 200 kpc (comoving) at SMG redshifts (z ~ 1 - 5). Thus as the massive, gas-rich, 
disks approach one another, both are contained in the beam for a significant time period. 
Since individually the disks form stars rapidly and are dusty enough to have high - but 
sub-SMG - sub-mm fluxes, two in the same ~ 200 kpc beam can yield a flux greater 
than the SMG threshold even if the disks' SFR is not elevated by gravitational torques. 
While this mechanism relies on merging so that the two disks are in the same beam, 
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Figure 1 . AzTEC 1 . 1 mm number counts predicted by our fiducial model includ- 
ing the contributions from the initial infall, first passage, and final coalescence of 
the two progenitor disks. We neglect the contribution from isolated disks as it is 
sub-dominant for the flux range plotted here. The ordinate is the cumulative num- 
ber of sources per square degree with 1.1 mm flux density greater than the value 
of the abscissa. The so lid line is the prediction calculated using the best-fit MF of 
Marchesini et al. (2009). The dotted lines show the range of predictions when the 
lcr random uncertainties and systematic uncertainties (from photometric redshifts 
and SED fitting) from Marchesini et al. are used. Uncertainties on our calculated 
fluxes (owing primarily to the sub-resolution ISM treatment and initial conditions 
of both the Gadget-2 simulations and the RT) are not included here but can have a 
significant effect on the predicted counts. See §H]for discussion. 

the sub-mm flux is not driven by the merger-induced starburst. As shown in Figure [Q 
in our model the combination of these two mechanisms can account for the observed 
number density of SMGs without invoking a non-standard IMF. 



4. Discussion 

One primary advantage our simulations is that they have the resolution necessary to 
model the gravitationally-induced tidal torques responsible for driving gas into the cen- 
ter of the merging galaxies and thus fueling the starburst and AGN. Furthermore, the 
high-resolution simulations enable us to perform RT on a realistic 3-D distribution of 
stars, AGN, and dust. We self-consistently calculate the thermal equilibrium tempera- 
ture of each dust species including dust self-absorption. These aspects of our models 
are essential for an accurate determination of the sub-mm flux. However, this accu- 
rate modeling is very computationally expensive, and we must model a large parameter 
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Figure 2. Same as Figure Q]but using only snapshots with gas fraction less than 
0.4. This choice eliminates most of the contribution from the infalling disks. 



space. The number of high-resolution models that must be run prevents us from per- 
forming cosmological simulations. One significant limitation that results is that ou r 
simulations do not include cosmological gas accretion (e.g.. iKeres et al.ll2005L 12009). 
Instead, we initialize the disks with high gas fractions (80%) so that a sufficient gas 
reservoir is present at the coalescence of the two galaxies. Our somewhat crude treat- 
ment of the galaxies' gas content likely causes us to underestimate the sub-mm flux at 
merger coalescence because the typical gas fractions at that time ( 15 - 30%) are less 
than the typical gas fractions observed (~ 40%) for both SMGs dTacconi et all 1 2006) 
and high-z massive star-forming galaxies of comparable mass (ITacconi et al.l 120100 . 
Higher gas fraction can lead to a stronger starburst, greater dust mass, and, conse- 
quently, higher sub-mm flux. On the other hand, inclusion of disks with gas fractions 
higher than those observed can bias the counts upwards. Thus our fiducial model may 
overpredict the contribution from the earliest stages of the merger when the two rapidly 
star-forming disks are both in the beam but not yet strongly interacting. 

To avoid including disks with unrealistically high gas fractions we have discarded 
all snapshots with gas fraction greater than 60% when calculating the counts plotted 
in Figure Q] To check the sensitivity of our prediction to our gas fraction cut, we also 
calculate the number counts using only snapshots with gas fraction less than 0.4. The 
results are shown in Figure [2] The counts are less than observed by a factor of a few 
for most bins when the best-fit Marchesini MF is used to compute the merger rates, but 
they are within the lcr upper limit except for the highest flux bin. 

Despite our relatively sophisticated treatment of the RT, significant uncertainties in 
the sub-mm flux calculation r emain. One of the most important is t he treatment of the 
sub-resolution ISM structure (Youn ger et al.1120091 : IWuyts et al.ll2009l) . As described in 
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$2j we assume that the dust density is smooth on scales below the Gadget-2 resolution. 
We have used the total gas/metal density in the Gadget-2 simulations to calculate the 
dust density, ignoring the SH03 multiphase breakdown. An alternate treatment of the 
sub-resolution ISM is to discard the cold phase of the SH03 model, effectively assuming 
the cold clouds have a negligible volume filling factor. Observations suggest that in 
local ULIR Gs the dense molecula r gas component of the ISM has a high volume filling 
factor (e.g.. ISakamoto et"ai1ll999h . leading us to adopt the former treatment, but this 
choice remains a source of uncertainty. To constrain this uncertainty, we have calculated 
the fluxes for a restricted set of models where we have discarded the cold phase, which 
typically accounts for > 90% of the dust mass. Predictably, this leads to hotter SEDs, 
which results in lower sub-mm fluxes (typically by a factor of a few). This provides 
unrealistically low attenuation when the attenuation from birth clouds is not included, 
but it is useful as a lower limit on the sub-mm flux from a simulated galaxy. Work to 
improve the Sunrise treatment of the sub-resolution ISM structure is underway. 



5. Conclusion 

We have, for the first time, combined high-resolution 3-D hydrodynamic simulations 
of isolated and merging galaxies, radiative transfer, and a semi-empirical merger rate 
model to predict number densities of SMGs. We find that isolated disks do not con- 
tribute significantly to the classical SMG population, whereas mergers contribute both 
because of the coalescence-induced starburst and because the two sub-mm bright but 
sub-SMG disks would not be resolved by the sub-mm instruments (e.g., SCUBA and 
AzTEC) used to determine wide-field number counts. We constrain free parameters by 
observations whenever possible and carefully test the effects of remaining uncertainties. 
We stress that we have not used a top-heavy IMF. Our predicted counts agree reasonably 
well with observations, though significant modeling uncertainties remain. The agree- 
ment is due to a combination of the relatively long sub-mm duty cycles of our simulated 
mergers (as much as a few hundred Myr for the mos t massive mergers, w hich is consis- 



ve merge 

tent with duty cycles estimated from observations; iTacconi et alj|2006l) . the relatively 
higher number density of massive, gas-rich mergers in our se mi-empirical model, and 
the reduction in counts observed by recent deep/wide surveys ( Austermann et al.ll2010l) 



relative to previous results. Our ability to match the observed counts without invoking 
a flat IMF indicates that, contrary to previous claims, observed SMG number densities 
do not provide sufficient evidence to reject the null hypothesis of a universal IMF. 
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